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The application of logarithmic transformations to speckle images is sometimes desirable in converting the speckle 
noise distribution into an additive, constant-variance noise distribution. The optical transmission properties of 
some bacteriorhodopsin films are well suited to implement such a transformation optically in a parallel fashion. 
I present experimental results of the optical conversion of a speckle image into a transformed image with signal- 
independent noise statistics, using the real-time photochromic properties of bacteriorhodopsin. The original and 
transformed noise statistics are confirmed by histogram analysis. 


Logarithmic image transformations may be useful 
any time the noise on an image can be described 
as multiplicative in nature. 1 This nonlinear oper- 
ation on the image separates the noise from the 
underlying object function by producing a new im- 
age with additive noise. Traditional linear filtering 
techniques can then be applied for purposes of im- 
age restoration, etc. One commonly occurring type 
of multiplicative image noise is speckle, found during 
imaging with spatially coherent light or radiation 
when the object of interest has a random surface 
roughness of the order of a wavelength and the imag- 
ing system cannot resolve the microscale of the ob- 
ject’s roughness. 2 In addition to allowing the easy 
application of linear filters for image processing, use 
of a logarithmic function has been recently shown 
to significantly improve object detection and recog- 
nition of speckle images through correlation with a 
matched filter after preprocessing the image with a 
logarithmic transformation to make the noise addi- 
tive and signal independent. 3 In this Letter I pro- 
pose and present experimental results of the optical 
logarithmic transformation of speckle images by us- 
ing the real-time photochromic properties of bacte- 
riohodopsin (BR). Although digital implementation 
is relatively easy, there may be situations such as 
the input of speckle images into an optical correla- 
tor in which it would be advantageous to perform all 
processing operations optically. BR films can accom- 
plish this transformation in real time and in parallel 
and are thus an attractive option. 

The statistical properties of speckle noise have been 
thoroughly studied, and it is well known that the im- 
age intensity generally follows a negative exponential 
probability-density function (pdf). 4 An image model 
that has been widely used to describe speckle images 
is the multiplicative model, 

r sp (x,y) = Ks(x,y)n sp (x,y), (1) 

where r sp (x,y) is the noisy speckle image intensity, 
s(x,y) is the incoherent intensity image of the ob- 
ject, and n sp (jt,.y) is the noise function. 15 We note 
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that Eq. (1) is strictly true only for spatially uniform 
areas of the object. 5 Assuming the model of Eq. (1) 
and negative exponential statistics for rc sp (x,y), it is 
easy to show that both the mean and the standard 
deviation of r sp (x,y) are proportional to s(x,y), which 
clearly makes the noise signal-dependent in nature. 
Because the logarithm of a product is mathematically 
equal to the sum of the individual logs, it is obvi- 
ous that a logarithmic transformation of an image 
described by Eq. (1) will change the noise from mul- 
tiplicative to additive. However, for speckle noise 
it also makes the new noise signal independent in 
nature. 6 This considerably simplifies the image pro- 
cessing for restoration purposes and also produces an 
image better suited for correlation with a matched 
filter for pattern detection. 

For digital image processing applications, the log- 
arithmic operation would be performed digitally. 
However, for some optical processing applications, 
such as correlation, it may be more efficient to also 
preprocess the images in this logarithmic fashion 
optically. We then require an optical material that 
exhibits a logarithmic transfer function. BR is an 
organically derived material that demonstrates many 
desirable properties for use as a real-time optically 
addressed spatial light modulator for optical pro- 
cessing applications. 7,8 It is a photochromic protein 
found in the photosynthetic system of a salt-marsh 
bacterium called Halobacterium salinarium. The 
BR molecule is contained in a membrane commonly 
called the purple membrane. The absorption of light 
initiates a photocycle in the BR molecule, which 
makes it a useful material as an optically addressed 
spatial light modulator. In the dark, the molecule is 
initially in the bR state, with an absorption spectrum 
peaked at —570 nm. On absorption of a photon in 
this band, the molecule quickly passes through sev- 
eral intermediate states by thermal relaxation until 
it reaches the M state, with its absorption spectrum 
shifted —160 nm towards the blue. Thus images 
may be written in a BR film based on the conversion 
of molecules from the bR state to the M state. In 
its native form, at nominal pH, relative humidity, 
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Write beam power density (mW/cm 2 ) 

Fig. 1. Transmittance of read beam as a function of write 
beam power density for wild-type BR film. Both beams 
have wavelengths of 514.5 nm. 

and ambient temperature, the lifetime of the M state 
is —10 ms. However, the image lifetime can be in- 
creased by orders of magnitude by the use chemical 8 
and mutagenic 7 methods. 

The film that we used in this experiment was 
prepared from the wild-type form of BR by Bend 
Research, Inc. The purple membrane material was 
extracted from the bacteria, suspended in a polymer 
matrix, and deposited upon an optical flat. The film 
was dried under controlled humidity conditions to 
lower the water content. The M-state lifetime was 
lengthened to — 10 s by means of chemical additives, 
producing a solution pH of approximately 9.0. Fi- 
nally, the film was sealed between the substrate and 
another glass plate. The film thickness was approx- 
imately 100 /xm, and the optical density of the film 
in the bR state at the peak wavelength was 2.4. 

The light transmission characteristics of the film 
were measured experimentally by a variable-strength 
write beam and weak read beam whose intensity 
was small enough so as not to affect the transmis- 
sion measurement. Both write and read beams were 
plane waves from an argon-ion laser at 514.5 nm. 
The film was first exposed to the write beam, which 
caused part of the bR state population to traverse the 
photocycle and switch predominantly to the M state. 
In a steady-state condition, the read beam was passed 
through the exposed area of the film and the trans- 
mitted intensity was measured with a photodetector. 
The data describing the relationship between trans- 
mitted intensity of the read beam, 7 read , and the inci- 
dent power density of the write beam, £ writc , are given 
in Fig. 1. It is clear that the function is logarithmic 
over the central part of the curve, i.e., 

ft ran s, read k lOgl-S^rj.jte) “ 1 “ C* , (2) 

for some constants k and e. Furthermore, the dy- 
namic range of the write power over which the trans- 
mittance of the BR film is logarithmic is large (>2 
orders of magnitude), which is required for the trans- 
formation of images with multiplicative noise such as 
speckle. Other films made from the wild-type form 
of BR also displayed behavior very similar to that 
shown in Fig. 1. 

I then performed an experiment to optically imple- 
ment a logarithmic image transformation with the 


BR film. The optical setup is shown in Fig. 2. The 
argon laser beam is collimated and then split into 
two beams, each controlled with an electronic shutter. 
With shutter A open and shutter B closed, a trans- 
parency in contact with a diffuser is illuminated and 
an imaging system forms the speckle image of the 
transparency on the BR film. The power level of the 
exposing beam is set to ensure that we are operating 
in the central logarithmic region of the transmission 
curve of Fig. 1. Immediately after the speckle image 
is recorded on the film, shutter A is closed, shutter B 
is opened, and the beam transmitted through the BR 
film is imaged into a CCD camera. The read beam 
power level is weak to ensure that it does not affect 
the image written in the BR film and also so that it is 
within the linear operating range of the CCD camera. 

The performance of the BR film was evaluated 
with a two-tone object encoded on the transparency. 
The right half of the object is clear, and the left 
half is a gray level. The speckle image of this ob- 
ject as obtained at the BR film plane is shown in 
Fig. 3(a). The image obtained by the CCD camera 
on transmission of the read beam through the BR 
film after recording is shown in Fig. 3(b). I analyzed 
both images in Fig. 3 and determined the first- and 
second-order statistics; these results are presented in 
Table 1. For the speckle image, the equality of the 
mean and the standard deviation in each half is con- 
sistent with the model of a negative exponential pdf. 
On the other hand, the image transmitted by the BR 



Fig. 2. Optical setup for logarithmic transformation of 
speckle image with BR film. Shutter A is opened dur- 
ing write process, and shutter B is opened during read 
process. BS’s, beam splitters. 



(a) (b) 


Fig. 3. (a) Original speckle image of a two-tone object, 

(b) logarithmically transformed image of the object ob- 
tained through BR film. 
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Table 1. Image Statistics (Absolute Numbers Are Unrelated between Two Images) 



Speckle Image 

Image Transmitted by BR 


Left Half 

Right Half 

Left Half 

Right Half 

Mean, p 

31.9 

98.5 

90.0 

122.0 

Standard Deviation, a 

29.0 

96.8 

37.4 

37.3 


Original speckle image 



Intensity value 

(a) 


8R transformed image 



Fig. 4. (a) Histogram of the right half of a speckle image, 

(b) Histogram of the right half of the image transmitted 
by BR film. The solid black curves represent the pre- 
dicted pdfs. 

film shows remarkably different statistics; the means 
of the two halves are different, but the standard de- 
viations are almost exactly equal. Thus the noise of 
the transmitted image is indeed signal independent. 
It can also be easily verified that the ratio of means of 
the two halves of the BR-transmitted image is consis- 
tent with that expected from the logarithmic trans- 
formation of the speckle image. 

Histograms of the speckle image and the BR- 
transformed image are presented in Fig. 4. The 
speckle image histogram closely resembles the ex- 
pected negative exponential pdf of our speckle model. 
The transformed image histogram clearly shows that 
the new noise distribution is additive in nature. For 
a logarithmic transformation as given in Eq. (2), the 


expected pdf of the transformed image noise n is of 
the form 

p N (n) = a I0 n/k exp(- y 10'"*) (3) 

for some constants k , a and y. The pdf predicted 
by Eq. (3) is slightly more asymmetric, with a longer 
tail to the left-hand side than the histogram data in 
Fig. 4(b), but the experimental results are neverthe- 
less quite similar to those predicted. The relatively 
minor differences can probably be explained by exper- 
imental noise and measurement error that tend to 
make the real distribution more Gaussian than the 
purely theoretical model. 

I have experimentally demonstrated the optical log- 
arithmic transformation of an image by using bacte- 
riorhodopsin film. In particular, the multiplicative, 
signal-dependent noise of a speckle image has been 
transformed into additive, signal-independent noise. 
Given that BR films have high spatial resolution, 
require no external processing, and may be cycled 
through the write/read/erase process millions of times 
without degradation, such films have the potential to 
function as a rapid, parallel optical logarithmic trans- 
former for images with multiplicative noise such as 
may be required in the operation of an optical corre- 
lator or in general preparation for subsequent optical 
image processing algorithms. 
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02-05-06. The author thanks Dan Smithey of Bend 
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